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Abstract — In this paper, we consider a robust lattice alignment 
design for K-usex quasi-static MIMO interference channels with 
imperfect channel knowledge. With random Gaussian inputs, 
the conventional interference alignment (lA) method has the 
feasibility problem when the channel is quasi-static. On the other 
hand, structured lattices can create structured interference as 
opposed to the random interference caused by random Gaussian 
symbols. The structured interference space can be exploited to 
transmit the desired signals over the gaps. However, the existing 
alignment methods on the lattice codes for quasi-static channels 
either require infinite SNR or symmetric interference channel 
coefficients. Furthermore, perfect channel state information (CSI) 
is required for these alignment methods, which is difficult to 
achieve in practice. In this paper, we propose a robust lattice 
alignment method for quasi-static MIMO interference channels 
with imperfect CSI at all SNR regimes, and a two-stage decoding 
algorithm to decode the desired signal from the structured 
interference space. We derive the achievable data rate based on 
the proposed robust lattice alignment method, where the design of 
the precoders, decorrelators, scaling coefficients and interference 
quantization coefficients is jointly formulated as a mixed integer 
and continuous optimization problem. The effect of imperfect 
CSI is also accommodated in the optimization formulation, and 
hence the derived solution is robust to imperfect CSI. We also 
design a low complex iterative optimization algorithm for our 
robust lattice alignment method by using the existing iterative 
lA algorithm that was designed for the conventional lA method. 
Numerical results verify the advantages of the proposed robust 
lattice alignment method compared with the TDMA, two-stage 
ML decoding, generalized Han-Kobayashi (HK), distributive I A 
and conventional lA methods in the literature. 

Index Terms — lattice codes, interference alignment, interfer- 
ence channel, MIMO, imperfect CSI 



L Introduction 

Interference is a fundamental bottleneck in wireless systems. 
This is partially due to the lack of understanding interference 
from an information theoretic view. For example, the capacity 
region of the two-user Gaussian interference channels has 
been an open problem for over 30 years. It was only shown 
recently that the Han-Kobayashi (HK) achievable region (Tl 
achieves the capacity region to within one bit |2 |. Lately, there 
have been some breakthroughs on the understanding of K- 
user interference channels. In |3] and the authors propose 
an interference alignment (lA) method to align interference 
onto a lower dimensional subspace of each receiver so that 



the desired signal can be transmitted on the interference- 
free dimensions. The authors show that lA is optimal in 
the degrees-of-freedom (DoF) sense (which is a high SNR 
performance measure) and the total capacity of the /-iT-user 
interference channels is given by f log(SNR) + o(log(SNR)) 
0. There are a number of extensions (|5], (|6] that have 
studied the application of lA in if-user quasi-static MIMO 
interference channels. However, the conventional lA method 
for iiT-user interference channels requires infinite dimension 
in time-varying or frequency-selective channels, and has the 
feasibility problem. For example, it is shown in [7] that 
conventional lA in quasi-static MIMO (A/ transmit and N 
receive antennas) interference channels is not able to achieve 
a per user DoF greater than . As a result, there is 

no satisfactory solution for quasi-static MIMO interference 
channels for large K due to the feasibility problem. 

Besides the feasibility problem, conventional lA solutions 
have assumed Gaussian input symbols and it is unclear 
whether it is optimal to employ Gaussian input for interference 
channels. With the Gaussian input symbols, the interference 
space is random and all the conventional lA methods try to 
align all interference to a smaller dimension space and utilize 
the remaining interference-free dimensions to transmit the de- 
sired signals. While Gaussian inputs (using random codebook 
argument) are capacity optimal in multi-access and broadcast 
channels JD, |l9l, they are not optimal for interference channels 
El, ifTOl . It is revealed in II 101 that structured interference is 
more preferred in interference networks. Instead of giving up 
the space that is reserved for the interference and transmitting 
the desired signal on the remaining interference-free dimen- 
sions, one could also exploit the structured interference space 
and transmit desired signals over the gaps. This motivates the 
lattices based study of interference alignment. In ifTTl . llT2l . 
the authors propose a lattice interference alignment method to 
align the interfering lattices on a common basis but it only 
works for symmetric SISO interference channels (where all 
cross links have the same fading coefficients) or a specialized 
class of 3-user SISO interference channels (where the products 
of the fading coefficients are assumed to be rational). In 
fT3l, fT4'|, the authors propose interference alignment over 
real line for X-user quasi-static MIMO interference channels 
with real channel coefficients and demonstrate that the DoF 
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— jj^fjK can be achieved. This approach is further extended 
to the complex channel in flSl for the compound MIMO 
broadcast channel. In addition to interference channels, lattice 
codes have also been widely studied in many communication 
networks [TSl such as point-to-point channels and wireless 
relay networks. Specifically, in [17], the author shows that 
lattice codes can achieve the capacity of ^ log(l + SNR) in the 
point-to-point additive white Gaussian noise (AWGN) channel. 
In |fT8|. the authors exploit the lattice codes and propose 
a compute-and-forward method in wireless relay networks. 
They show that each relay can decode a linear function of 
messages from multiple source nodes by using lattice codes. 
A destination node can decode the desired message, given 
sufficient linear combinations forwarded from the relay nodes. 

All these results indicate that we can potentially benefit from 
the structured interference created by lattice codes. However, 
in order to have a practical method for exploiting the advantage 
of the structured interference due to lattice transmissions, there 
are still some key technical challenges to be addressed. 

• How to align the received lattices on a common basis 
over irrational quasi-static MIMO interference chan- 
nels: Creating an ideal structured interference space 1, 
lfT2l requires that all the interfering lattices to be received 
on a common basis at each of the K receivers. However, 
this is a difficult requkement and is impossible for 
irrational fading matrices. While JTT], ifTZl study lattice 
alignment methods for K user quasi-static interference 
channels, the methods therein only work for symmetric 
interference channels or a specialized class of 3-user 
SISO interference channels. These approaches cannot be 
used for general non-symmetric in^ational interference 
channels. 

• How to exploit structured interference at finite SNR: 
In |fT3l , ||T4|, the authors propose a lattice alignment 
algorithm for MIMO interference channels, which is 
optimal in the DoF sense. However, it is not known 
whether this approach can be applied at finite SNR, which 
is an important operating regime in practice. In llT9l . the 
authors propose an ergodic interference alignment method 
which could work at finite SNR regime. However, it 
requires either infinite time or frequency extension and 
cannot be applied for constant channels. 

• How to ensure robustness due to imperfect knowledge 
of CSI: All of the above works assume perfect knowl- 
edge of channel state information (CSI) for facilitating 
interference alignment. In practice, this is not possible 
and interference alignment performs very poorly with 
imperfect CSI |20|. It is quite challenging to incorporate 
robustness with respect to (w.rt.) imperfect CSI and 
exploit the structured interference space at the same time. 

In this paper, we explore a robust precoder and decorre- 
lator design for A'-user general irrational quasi-static MIMO 
interference channels with imperfect CSI. We propose a robust 
lattice alignment method (which does not have the feasibility 
problem) for general irrational constant interference channels, 
and deduce a two-stage decoding algorithm to decode the 
desired signal from the structured interference space. We 



derive the achievable data rate based on the proposed method 
and formulate the precoders, decorrelators, scaling coefficients 
as well as the interference quantization coefficiently] design as 
a mixed integer and continuous optimization problem 11211 . By 
utilizing the alternating optimization technique and by exploit- 
ing the structured interference, we derive a low complexity 
algorithm to determine the precoders, decorrelators as well as 
the interference quantization coefficients. The effect of imper- 
fect CSI is also accommodated in the optimization formulation 
and hence, the derived solution is robust to imperfect CSI. 
To illustrate the benefit of the proposed method, we compare 
the achievable data rates with those of the conventional lA 
method lO, the distributed I A method based on alternating 
optimization fSl, the TDMA method, the brute-force two- 
stage maximum likelihood (ML) decoding method and the 
generalized HK method jll, ||2| with random Gaussian input 
symbols. We show that the proposed method achieves the same 
DoF as the conventional lA method if the problem is feasible. 
On the other hand, when the conventional lA is not feasible, 
our proposed solution can still offer significant performance 
gain compared with these baselines. 

This paper is organized as follows. In Section |ll] we 
outline the system model of the ii'-user quasi-static MIMO 
interference channels. In Section Hill we review some common 
interference alignment methods in the Uterature and provide 
some preliminary discussions on the lattices. In Section |IV] 
we propose the robust lattice alignment method and formulate 
it as a mixed integer and continuous optimization problem. In 
Section [V] by using the alternating optimization technique, 
we derive a low complexity solution. In Section IVII we 
derive closed-form analysis of the proposed method and the 
baseline methods under specialized channel realizations. The 
performance of the proposed method and the existing methods 
under complex channel realizations is illustrated in Section 
IVIII We conclude with a brief summary of the results in 
Section IVm] 

II. System Model 
A. K-user Quasi-Static MIMO Interference Channels 

We consider the ii'-user quasi-static MIMO Gaussian in- 
terference channels as illustrated in Fig. [T] Each transmit- 
ter, equipped with M antennas, tries to communicate to its 
corresponding receiver, which is equipped with N antennas. 
Specifically, we consider a block of T channel symbols. The 
channel output at the fc-th receiver is described as follows: 

Yfe =^Hfc,X, + Zfc, (1) 

i 

where, 11^:^ is the N x M MIMO complex fading coefficients 
from the i-th transmitter to the fc-th receiver Xi is the M xT 
complex signal vector transmitted by transmitter i. Zk is the 
N X T circularly symmetric AWGN vector at receiver k. We 
assume all noise terms are i.i.d. zero mean complex Gaussian 
with unit variance. 

' Please refer to Section IIVI for details about these design parameters. 
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Fig. 1. Quasi-static complex X-user MIMO Gaussian interference channels. 
Each transmitter, equipped with M antennas, tries to transmit L independent 
data streams to its corresponding receiver, which is equipped with N antennas. 



B. Precoding at Transmitters 

In this paper, we assume that each transmitter transmits 
L independent data streams. Specifically, v' is the M x 1 
precoder for the l-th complex data stream x' (1 x T) at 
transmitter i. Therefore, Xi = 'l2i'^Wi- The average power 
for each data stream is given by yE[||x'|p] = P. The 
precoders are chosen such that the total transmit power for 
each transmitter is no more than jP, i.e.. 



El 



where llv'lP 



J\\2 



(2) 



C. Imperfect Channel State Information Model 
The imperfect CSI model can be expressed as: 



(3) 



where H^i is the estimated CSI that is known at the central 
coordinator. A^^ is the CSI error satisfying: 



{A. 



Ai 



{Afe, :Tr{AH(Afc,)^}<e'}. (4) 



This deterministic channel estimation error model is widely 
used in the literature ['2T|-|'24'|. Therefore, given the estimated 
CSI Hfci, the uncertainty of Hfei is an ellipsoid centered at Hfci 
with radius e. 

III. Summary of the Existing Interference 
Alignment Methods and Lattices 

A. Interference Alignment on a Quasi-Static MIMO Space 

It is shown in that in the 3-user quasi-static MIMO 
interference channels with A/ > 1 antennas, ^ DoF can 
be achieved. Specifically, when M is even, each transmitter 
sends ^■ independent streams over the M x precoder 
for transmitter i, i.e.. 



M/2 



(5) 



Hi2 ^ H23 ^ Hzi 

H21 H32 Hi3 



(6) 



The precoders {Vi}^^^ are designed in terms of the channel 
matrices, where the dimension of the interference space is 
equal to M/2 at each receiver. Each receiver can simply 
cancel the interference by zero-forcing and then decode the 
desired M/2 streams. When M is odd, similar results can be 
obtained by considering a two time-slot symbol extension of 
the channel, with the same channel coefficients over the two 
symbols. However, it is shown in [7| that the above method 
is not able to achieve a per user DoF greater than ^^"^^ . 

B. Interference Alignment on a Lattice 

In [T2I, the authors consider lA on a lattice for some 
special 3-user SISO interference channels {M — N = I). 
Specifically, the real channel fading coefficients have to satisfy 
the following requirement: 

p 

q 

where p and q are integers such that gcd{p,q) = 1. gcd{-) 
is the greatest common divider of the integers. Transmitter 
i chooses the lattice Aj generated by using construction A 
described in ||251 . In order to align the interfering lattices at 
each receiver, the lattices are chosen by: 

H12A2 = PH13A3, H21A1 — qH23A3, H31A1 = 11^,2^2- 

(7) 

As a result, a two-stage decoding algorithm can be done 
at each receiver, which decodes the aggregate interference 
first and then subtracts the aggregate interference from the 
received signal to decode the desired information. Note that 
the symmetric channel considered in [11] is a special case 
in (|6]l. This method requires strong interference channels to 
decode the interference first (treating the desired signal as 
additive noise). Furthermore, the requirement of adding the 
interfering lattice on a common basis at each receiver as in 
(|7]i becomes infeasible for general irrational X-user quasi- 
static interference channels. A simple counter example is given 
below. 

Example 1 (A Simple Counter Example): A simple 
counter example is given by: yi — xi + \/2x2 + 
a/3x3 + zi,y2 = v^xi + X2 + V7x3 + Z2,y3 = 
\/lTxi + vT3x2 + X3 + Z3. It is not possible to align 
the interfering lattices on a common basis at each of the 3 
receivers. ■ 

C. Interference Alignment over the Real Line 



It is shown in |fT4| that 



MN 
M+N 



K DoF can be achieved 
for quasi-static MIMO using lA on the real line lfT3l . This 
approach is extended to the complex channel in flSl for the 
compound MIMO broadcast channel. Here we review the basic 
idea for the SISO interference channels in llT3l . and it is easy 
to extend to the MIMO case as in |14l, (l3. Specifically, the 
received signal at receiver A; can be represented as 1.13] 



yk = A\ y \ HkkTkiXk 



1=0 



^ 1 

E E. 

i—l^i^k 



HkiTiix 



Zk, 



(8) 
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where A controls the input power of all users, x'j, G {~QtQ)z 
is one of the integers in the set (— Q, Q), and carries infor- 
mation for the Z-th data stream for user fc. Tki is a constant 
real number, which is the direction of the transmitted ^-th data 
stream. It is chosen as monomials with variables from channel 
coefficients. Ik is the aggregated interference caused by all 
users. The data streams are aligned if they arrive at the same 
direction, e.g., interference x\ and X2 are aligned at receiver 
3 if i/aiT^ = H:i2Ti. By carefully desi ening the transmit 



MN 
M+N 



K DoF can be 



directions, it is shown in ifTJl . lfT4l that 
achieved. However, this method requires infinite SNR and it 
is not known if this method could be modified for finite SNR 
regime. 



D. Lattices 

In this section, we shall review some preliminaries on 
lattices from ||T6l , ifTSl and the references therein. A T- 
dimensional lattice A is a set of points in R^, and is given in 
terms of the lattice generator matrix L e R-^^-^: 



A = {x = Lw : w e Z^}. 



(9) 



A lattice quantization is a function, Qa : R-^ — > A, 
that maps a point x to the nearest lattice point in Euclidean 
distance: 



(5a(x) = argmin||x- A||. 

AeA 



(10) 



The fundamental Voronoi region, V, of A, is the set of points 
in R-^ closest to the zero vector, i.e., V = {x : Qa(x) = 0}. 
The modulo-A operation w.rt. the lattice is defined as 



X Mod A = X - Qa(x), 



(11) 



which is also the quantization error of x w.rt. A. 

A nested lattice code Z\ is the set of all points of a fine 
lattice Ai that are within the fundamental Voronoi region V 
of a coarse lattice A: 



£1 = Ai n V = {x : X = A Mod A, A e Ai}, 



(12) 



where A is said to be nested in Ai, i.e., A C Ai. Please refer 
to Fig. 2 in [18J for an illustration of nested lattice. The rate 
of a nested lattice code is given by: 



«.ilog(|£.|, = i,„g^, 



(13) 



where Vol(V) is the volume of V. 



IV. Robust Lattice Alignment 

In this section, we introduce the framework of our ro- 
bust lattice alignment method and the associated two-stage 
decoding algorithm. We derive the minimum achievable data 
rate and design the precoder, decorrelator and the interference 
quantization coefficients via an optimization approach. 



A. Encoding at Transmitters 

In this paper, we adopt the lattice encoding method, which 
is used both in point-to-point channels 1 17| and relay networks 
iiSJ. Specifically, the data stream x'j. is given by: 



Mod A + j 



Mod A, 



(14) 



where A G R^ is the coarse lattice used in the nested 
lattice encoding method for the transmitted symbols xj,, Vfc, /. 
{t'j., t^} are points in the nested lattice code L\ (corresponding 
to the fine lattice A'^) that carry information, {dj., dj,} are the 
dither vectors fTTl , |[T81 which are independently uniformly 
distributed over V and are available to all transmitters and 
receivers. Specifically, all the nested lattices {Aj,} are AWGN 
good, and the coarse lattice A is quantization goocd In ifTSll . 
such lattices and dither vectors satisfying the transmit power 
constraint yE[||x^|p] = P are proved to exist. 

Remark 1 ( Complexity of the Encoding Method): The 
complexity of the encoding method is similar to that in the 
baselines. In the baseline methods, standard Gaussian random 
codebook is assumed to achieve the mutual information rate. 
On the other hand, there exists efficient coding schemes 
by using a scalar constellation coupled with a linear code 
and this could come close to the nested lattice performance 
(without the shaping gain) ll26ll . ■ 

Note that, dithering is a common randomization technique 
in lattice quantization for source coding |17|. The following 
two lemmas from iflTl . ifTsl capture the key properties of the 
dithered nested lattice codes. 

Lemma 1 (Erez-Zamir HI 71/ ).- For any random variable t e 
V, if d is statistically independent of t and uniformly dis- 
tributed over V, [t — d] Mod A is uniformly distributed over 
V, and is statistically independent of t. ■ 

Lemma 2 (Nazer-Gastpar SWU): Let z ~ A/'(0,Itxt) and 
d" be statistically independently uniformly distributed over 
V with P = ^E[||d;^||2], and = az + E».„ ^-"df for 
some constant a and 0". The density of zj. is upper bounded 
by the density of an i.i.d. zero-mean Gaussian vector z'j, 
whose variance approaches Nl as T — 00, where Nl — 

Remark 2 (Interpretation of the Lemmas): Lemma [T] as- 
sures that the input power exactly meets the power constraint 
iflTl . Lemma |2] indicates the non-Gaussian noise z'j. is nearly 
Gaussian as the code length T increases. ■ 

B. Lattice Alignment with Imperfect CSI 

In this section, we shall discuss the proposed robust lattice 
alignment method (using vector space strategies) and the 
motivations. When Gaussian signals are transmitted in interfer- 
ence channels, the associated interference space is random as 
illustrated in Fig. |2]rT0|. Intuitively, the aggregate interference 
will fill the entire signal space and there is no room left 
for desired signal transmission. Furthermore, the penalty of 
interference scales with the number of interferers. As a result, 

-please refer to |16| for the definition of goodness for lattice codes. 
Specifically, the lattice codes should have both good statistical and good 
algebraic properties 1181 . 
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Fig. 2. Illustration of random interference (due to random codes) and 
structured interference (due to structured codes) in interference channels 
1101 Fig. 4, Fig. 5]. The resulting random interference covers the entire 
space, preventing the desired receiver from decoding. On the other hand, the 
structured interference allows the desired receiver to decode between gaps. 



the key idea behind the lA methods in lH, H is to aUgn all 
interferences to a lower dimensional subspace and utilize the 
remaining interference free dimensions to transmit the desired 
signals. However, one problem associated with the lA method 
in lO, m is the feasibility problem. In fact, it is shown in fl] 
that the per user DoF greater than is not achievable. 

On the other hand, when structured lattices are used at 
the transmitters, the aggregate interference at each of the 
receiver may be structured (if the interfering lattices are 
properly aligned) as illustrated in Fig. l2l llT0l . In this case, the 
interference space contains regular gaps that can be utilized 
to transmit desired signals. More importantly, the penalty of 
the interference does not scale with the number of interferers. 
Therefore, another possible direction of alignment is to align 
the transmit lattices on a common basis at each of the K 
receivers to create a structured interference space as illustrated 
in Fig. [2] |10|. However, it is not always possible to simulta- 
neously align all the interfering lattices, as shown in Example 
[T] Furthermore, perfect lattice alignment requires perfect CSI, 
which is impractical. 

Motivated by the advantages of having structured interfer- 
ence, we propose a robust lattice alignment method in which 
the precoders are designed to align the received lattices as 
much as possible. We accept the fact that lattice alignment may 
not be perfect (due to infeasible irrational channel coefficients 
and imperfect CSI) and we try to minimize the effect of 
the residual lattice alignment error. Specifically, instead of 
decoding and subtracting each interference individually, we 
wish to decode the structured aggregate interference composed 
of all the interferences and remove them all at once. Further- 
more, due to the structural constraints of the lattice, we must 
choose integer coefficient for each data stream to approximate 
the equivalent channel coefficient. Therefore, the aggregate 



interference lattice for the data stream is written as: 



(15) 



where a" G Z + jZ (complex integer) is defined as the inter- 
ference quantization coefficient to approximate the equivalent 
channel coefficient for x", and a'^, = (see Fig. |3] for an 
illustration). Specifically, given the precoder v" for x", the 
decorrelator at receiver fc, and the interference quantization 
coefficients {a"}, the alignment error is given by: 



(16) 
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Fig. 3. Illustration of residual alignment error after stage I decoding 
and how it is affected by different choices of stage I decorrelator uj. and 
interference quantization coefficients a = Specifically, a" is 

the interference quantization coefficient and Hjj^v" is the equivalent channel 
for the stream x" respectively. 




Fig. 4. The block diagram of the two-stage decoding algorithm for decoding 
the desired signal xj. . Specifically, is the observed signal for user k given 
by )19K ij, is the structured aggregate interference given by ( 1211 , uj. is the 
first stage decorrelator, aj. are the interference quantization coefficients, cj. 
is the scaling coefficient, and uL is the second stage decorrelator. 



Remark 3 (The Effect of Design Parameters on 1,,): Note 
that if the channel coefficients {ilki}fk=i irrational, the 
alignment error le may not be zero but the design parameters 
{ u[. , vf , af } in (fTSI l can be chosen to minimize the effect of 
the alignment error as illustrated in Fig. |3] We shall elaborate 
the precise optimization problem in Section IIV-EI ■ 

Based on the aggregate interference lattice in (115b . a two- 
stage decoding algorithm can be constructed at each of the 
receivers as illustrated in Fig. E] The first stage is to decode 
the structured aggregate interference, whereas the second stage 
is to decode the desired signal after canceling the decoded ag- 
gregate interference. The decoding process and the associated 
error analysis in the presence of the residual alignment errors 
are discussed in detail in the following subsections. 



C. Stage I - Decoding the Aggregate Interference 

Note that the received signal at the fc-th receiver is given 
by: 

Yfc = Hfcfcvlxi + E HfcfcV^x^ + ^ ^ Hfe^v^xr + Z,. 

(19) 
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log 



oo 



l|uLlP+^i: l«)"Hfc,v7-a7l+e||v,"|M|u'J| 



if a? ^ 
if a" = 



(17) 



4 < A^L = log (- 



P 



(18) 



At stage I, a complex linear decorrelator u\ (N x 1) is 
applied at receiver fc, i.e., 



desired signal 



+ ^(4)^H,,;vrxr+K)^z, 



inter-stieam interference 



inter-user interference 



(20) 

From y{ (IxT), we wish to decode the structured aggregate 
interference, i.e.. 



5R{a^xn modA+j y Sja^x^} mod A, 

(21) 

where e Z + jZ, and a{ - 0. Let = {{anfLjLi 
denote the interference quantization coefficients for decoding 
the desired data stream xj.. 3?{-} and denote the real 

and imaginary parts, respectively. To successfully decode l'^., 
there are some requirements on the transmit data rates of the 
streams in the term. The achievable rate region of the K 
users for successful stage I decoding under the imperfect CSI 
model in (|4|i is summarized in Lemma |3] 

Lemma 3: (A Sufficient Condition for Successful Stage I 
Decoding under Imperfect CSI): Under the imperfect CSI 
model in the structured aggregate interference I^. can be 
decoded from with arbitrarily small error probability if 
the data rate i?" satisfies the condition (for all i, n) given in 
([TtI i. where ||vf || • ||u[.|| means the product of the two terms 
{||v?||,||ui||}. ■ 
Proof: Please refer to Appendix |A] Note that a similar 
result with perfect CSI is derived in lITSl for single-stream 
wireless relay networks. 



Note that when 



~ for all {i,n}, i.e., 



there will be no stage I decoding for the desired data stream 
x^ and the proposed method reduces to the conventional 
precoder-decorrelator optimization with one-stage decoding. 
The first term and the second term in the denominator of (ITTi 
correspond to the noise term and the residual lattice alignment 
eiTors, respectively. 



D. Stage II - Decoding the Desired Signal 

After successfully decoding the structured aggregate inter- 
ference at stage I, we wish to decode the desired signal x^. 
at stage II. This is illustrated in Fig. |4] Specifically, the scaled 
structured aggregate interference c\\\ is subtracted from the 
received signal and the desired signal is decoded via y\ given 



by: 



yi 



u - (5^)^Hfc„,vlxl If a', = 5t 

SR{(at)«Yfc-4lt}]ModA+ 

4s>{(aL)"Y,-4iL}]ModA 



otherwise 



, (22) 



where is the second stage decorrelator, cj^, G Z + jZ is the 
scaling coefficient for the l-\h data stream at receiver fc, 5^ 
is the (5-function (all a" = except aj^ — 1), and 

T' }] Mod A 

Mod A 



-.1 \H 



(a) 



5R{(Si)^Y,}-5iK}5R{li}+3{ci}S{li} 
5R{(Si)«Y,.}-[KK}K{li}]ModA 



-[3{4}3{li}]ModA 



Mod A 



3«{E.„[(5i)^HHvf-4ar]xr 



Mod A 

(23) 

gi]Mod A) + 

g2]Mod A for all gi,g2 S M^, (b) follows from c[ e Z + jZ 
and [c([gi]Mod A)]Mod A = [cgijMod A for all c e Z. 
Similarly, we have 



where (a) follows from [gi + g2] Mod A 



3{(2l.)^Y,-4li}]ModA = 



Mod A. 



(24) 



From ( |22] | note that, if we set aj^. = Sf^^, then l[, 
and hence the stage II decoding in ( l22b would try to directly 
null off the interference x^. Similarly, the sufficient condition 
in terms of the rate region for successful stage II decoding is 
given in the next lemma: 

Lemma 4: (A Sufficient Condition for Successful Stage II 
Decoding under Imperfect CSI): A sufficient condition for 
successful stage II decoding under the imperfect CSI model 



_ Q in dUi is given in dTSl l. where is an indicator function. 



Proof: The proof follows from a similar approach as in 
Appendix |A] Specifically, it is obtained by replacing yj, and 
ijj (obtained from (|23] | and HtM ) with yj, and x'j, in the proof 
of Lemma [3] respectively. Note that, from (l23T l and (l24t . the 
equivalent channel coefficient for data stream x" in stage II 
decoding is (u'j.)-^Hj,jv" — c'j.a". Furthermore, since we wish 



to decode x' 



1 • X 



fe 



I 



1 



{i—kSz7i—l} 



indicates the coefficient for x), and the coefficients for the 
other data streams are all 0. ■ 
Note that the rate region in (fTsT l is a function of inter- 
ference quantization coefficients a, scaling coefficients c = 

^ We may need to scale the decoded aggregate interference first before 
subtracting from the received signal so as to compensate for tlie amplitude 
change due to the second stage decorrelator u^.. 
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{{CfelfeLilz^i' precoders v, and second stage decorrelators 

Remark 4 ( Complexity of the Decoding Method): The de- 
coding complexity of the proposed design and the baselines 
are similar For instance, all baseline methods have assumed 
ML decoding to achieve the mutual information rate. On the 
other hand, there exists efficient lattice decoding methods ifTTl . 
|[T8l that could exploit the lattice symmetry in the decoding 
process. ■ 

E. Precoder, Decorrelator and Interference Quantization Co- 
efficients Optimization 

In this subsection, we shall formulate the precoders, decor- 
relators, scaling coefficients and interference quantization co- 
efficients design problem formally as an optimization problem. 
The problem consists of the following components: 

• Optimization Variables: Specifically, the optimization 
variables consist of the set of precoders v, the set of 
stage I decorrelators u, the set of stage II decorrelators 
u, the set of interference quantization coefficients a (used 
in stage I processing), and the set of scaling coefficients 
c (used in stage II processing). 

• Optimization Objective: For fairness, we consider the 
worst-case data rate as the optimization objective, i.e., 

max i?min, where = miiifc^; (^^., 

u,u,v,a,c 

• Optimization Constraints: The optimization constraints 
consist of the transmit power constraint X^^i 11"^*: IP — 
7, Vfc, the interference quantization coefficients constraint 
ajj. e (Z+jZ)^^, and the scaling coefficients constraint 

4 e z + jZ. 

As a result, the optimization problem is summarized below: 

arg max m\n{^i{,']l{) 

u,u.v,a.c L,K 
L 



u ,u ,v ,a ,c } 



s.t 



1=1 



rl l|2 



< 7, Vfc 



aj. G (Z + jZ)^^; 4 e Z + jZ 
(25) 

The above optimization problem involves complex 
{u, u, v} and integer variables {a, c}, which is a mixed 
integer and continuous optimization problem. As a result, the 
problem is non-convex and requires exhaustive search [21]. 

V. Low Complexity Iterative Solution 

In this section, we shall propose a low complexity iterative 
algorithm by exploiting the special structure of the problem 
in 



A. Properties of the Optimal Interference Quantization Coef- 
ficients {a*} 

Although obtaining the optimal integer solution {a*} in (l25l l 
is a difficult problem, we have the following lemma to reduce 
the search space for a*. 

Lemma 5 (Properties of the Optimal a*}: The optimal in- 
teger solution a* in (IZST i should belong to the following set, 
i.e.. 



(ai)*e^ = {a^^ 
Vr e (Z + jZ) 



^Z + jZ,Vaf 7^0, 
and |r| 7^ 1} \/k,l. 



(26) 



Proof: Suppose {K)*, K)*, K)*, (c' )*, (a',)* ^ A} 



are the optimal solution. Suppose that 



e Z + jZ, 



and \r\ > 1. Although a'j. will influence /i'^ in ( flTl l and 
in ([18]), it is easy to verify that /"K^^^, ■^^^) > 

/.i(K)*,(ai)*), and ^^(4)*r, (41!) = /^^Vfe)^k.)*)■ 
In other words, a higher data rate is achievable with the 
designed parameters {(v[)*, (Qj^.)*, (4.)*r, -^j. As 

a result, {(v^)*, (H^)*, K)*, (4)*, (a^)* ^ A} cannot be 
an optimal solution. ■ 

B. Optimization of {u, u, c} under fixed {v, a} 

In this section, we fix the precoder v and interference 
quantization coefficients a, and we optimize the remaining 
parameters {u, u, c}. Note that, uj, only influences /ij, in ( [TT] ), 
and uj, only influences Jl^ in ( fTSl ). For given c\., we shall first 
determine the optimal \i\ and vl\ by maximizing /i^ and /Ij,, 
respectively, which can be obtained by solving the following 
convex problems 



~ argmin 

|lK)^H,,vf-af| 
(uj,)* = argmin ( ||u' 

|l(ai)^H..v 



I|2 



cia" 



l{i=fc&n=/}l + 



u . 



(27) 

For perfect CSI (e = 0, H = H), closed form solutions 
exist for (|27| i and the optimal solutions are given by li27l : 

-1 



1 
p 
1 

p 



(28) 



(29) 



where W - [Hfeiv}, H^ivf, Hfesv^ 



KL X iV matrix, aj, = [aj, af , 02, .. 

■1 Cfc*^! ' 



2i ■ 



••,HfcKv|^]^ is a 
is a KL x 

..,1 



c\aW^ IS 



1 vector, and /S),. = \c\a 
also a KL x 1 vector On the other hand, when the CSI is 
imperfect, there is no closed form solution. Since (l27t is a 
standard convex problem, (u^* and (i3^)* can be obtained 
iteratively using an interior-point method or efficient gradient 
search 



Next, we shall optimize the complex integers c'j, for a given 
u'j.. The solution of (4)* for given (u[.)* is summarized in 
the lemma below. 

Lemma 6 (Optimal Integer Solution of (cj,)* given (u^)*j.- 
Given (u'j.)* in ( |27] ). the optimal (4)* is given by: 



(4)* = argmin /(4), 

!R{4}e [sR{r},!R{K}] .Oic'js [s{t},S{k}] 

where 



(30) 



r = (4)* - (1 + j) and k = (4)* + (1 + j), (31) 
/(4) = 

E i(2i)^HHvr - 4«r - i{.=/c&„=/}i + eiivrii • 



(32) 



8 



and 

(21)* =argmin/(24). (33) 

Proof: c\ only influences pt^. in ( fTSl l. If we relax c\ E 
Z + jZ to c'j, G C, the function /(c*^) is a convex function. 
Therefore, if (c^)* = argmin^^f^gj- /(c*^), the optimal is 
one of the complex integers close to (c^.)*, i.e., 

{cir= argmin /(ci), (34) 

5R{c'Je [5?{r},SR{K}] ,aK}e [3{r},a{K}] 

where r = {c[r - (1 + j), « = + (1 + j)- ■ 
As a result, given {v, a}, we shall use the following 
algorithm to optimize {u. u, c}. 

Subalgorithm A: Optimization Algorithm for {u, u} and 
c under fixed {v, a} 

> Step 1: For a given realization of v, £ A, initialize 
c(0). Set the iteration steps m = 0. 

• Step 2: For the given {v, a}, solve the convex op- 
timization problem (|27] | to obtain the corresponding 
optimal first stage decorrelator u. 

• Step 3: For the given {v, a, c(to)}, solve the convex 
optimization problem (l27T i to obtain the correspond- 
ing optimal second stage decorrelator u(m + 1). 

• Step 4: For the given {v, a, u(to + 1)}, obtain the 
optimal positive integers c(m + 1) by solving the 
problem ( [30l l. Specifically, solve the convex problem 
( [33] l to obtain the relaxed value (c^)*, and then check 
the nearby integers around (c^)* to obtain c(m + 1) 
as shown in ( l30l l. 

• Step 5: Continue until u(m) = u(m + 1), and 
c(m) = c(m + 1). 



The convergence proof is given in Appendix |B] 

C. Optimization o/{v,a} under fixed {u, u, c} 

Given {u,u, c}, optimizing {v,a} is not trivial. To obtain 
low complexity solutions, we shall first relax the interference 
quantization coefficients from a^ e (Z+jZ)^^^ to a[. e C^^. 
Define 

5i(v,ai) = ||ui|p+FE,„ 
|K)^HHV?-a?| + e||v?||.||4|||' 
5i(v,ai) = ||Gil|2+F^^_„ 

\iu[)"n,,vf - c[af - l{,=fe&„=i}| +e||vr|| • ||i3i|||', 

where the right hand side of these two equations are the 
denominator of (fTTI l and ( fTSl l. respectively, {fffc, 5^.} are convex 
functions w.r.t. {v, aj^.} for all {l,k}. Given {u,u, c}, the 
max-min problem in (IZSl l is equivalent to 

maxv,a.t t 

s.t. t < A t < Vfc, / 

L 

ElK-ll'<7,Vfc, (35) 
1=1 



which is not a convex problem, because log (^[(v, a'^,)) and 
log (g[.(v, ajj,)) are not convex functions w.rt. {v, a^} for all 
{l,k}. However, we show that (|35] l is equivalent to a convex 
problem given in the following lemma. 

Lemma 7 (Equivalent Convex Problem): The problem in 
(l35] l is equivalent to the following convex problem 

minv,a,t t 

s.t. 5i(v,ai.) < t and5[(v,a^) < i,V/,fc 

L 

1=1 

Proof: We show that given {u, u, c}, the max-min prob- 
lem in ( |25] ) is equivalent to 

argminmax(5[,5[). (37) 

v,a 

First, note that given {v, a}, the optimal 
{P,fc*} obtained from argmaxi_fe(g[, 5^) is the 
same as argmini_fc(/i^, /I^). Furthermore, note that 
arg minv,a(5r* > dk* ) = arg maxv,a(/ifc. , mC- )■ 

Therefore, 

argminmax(5[,5[) = argmaxmin(^^, /I^). (38) 

v,a i,/c v,a i,fc 

Similarly, we can also verify that (|38] l is equivalent to (|36] l. 

■ 

As a result, given {u, u, c}, we shall use the following 
algorithm to optimize {v, a} by using an interior-point method 
|28, Chap.ll]. 



Subalgorithm B: Optimization Algorithm for {v, a} 
under fixed {u, u, c} 

• Step 1: For a given realization of {u, u,c}, initialize 
g > 0, {vq, ao}, V > \ and tolerance <; > 0. 

• Step 2: Solve the unconstrained convex problem 
{v*(q),a*((7)} = argmiuv^agt - Efcj(log(i - 
fffe(v, ^i)) + log(i - g[(v, a'fe))), starting at {vq, ao} by 
efficient gradient searchQ 

. Step 3: If i < stop and set {v*, a*} = {v*(g), a*(g)}, 

otherwise go to the next step. 
. Step 4: Set {vo,ao} = {v*(g), a*(g)}, 

• Step 5: Set q = vq, and go to the step 2. 



The convergence proof is shown in Appendix [B] 



D. A Summary of the Overall Solution 

The top-level optimization algorithm is summarized below 
(and illustrated in Fig. |5]i; 

For example, using the Newton's method f28' Alg 9.5]. 
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lniliali«ilion 



v(0),a(0),c(0) 



Optimize u by solving 
problem (27) 



Optmiize u by solving 
problem (27) 



11 (ra) 



Optimize c by solving 
problem (30) 



Subalgorithm A 



m = m + l 




Optimize v_ a by 
Subalgorithm B 



I v(H).a(n) 



n = n + l 




Fig. 5. Illustration of the top-level optimization algorithm. Specifically, v 
denotes the precoders, a denotes the interference quantization coefficients, u 
denotes the first stage decoiTelators, u denotes the second stage decorrelators, 
and c denotes the scaling coefficients. 



Algorithm 1: Top-Level Optimization Algorithm 

. Step 1: Initialize v^(0), a'fc(O) e A and 4(0). Set 
the iteration steps m — 0. 

• Step 2: For the given {v(to), a(TO)}, apply the subal- 
gorithm A to obtain {u(rn + 1), u(to-|- 1), c(m-|- 1)}. 

• Step 3: For the obtained {u(m + l), u(m-l-l), c(m-|- 
1)}, apply subalgorithm B to obtain the correspond- 
ing optimal {v{m + 1), a(m + 1)}. 

• Step 4: Continue until u(m) = u(m + 1), u(m) = 
u(m + l), c(m) — c(m + l), v(m) ~ v(m + l), and 
a(m) = a(m + 1). 

The convergence proof of Algorithm 1 is shown in Ap- 
pendix |B] Finally, suppose (aj,)* is the solution to Algorithm 
1, we could quantize (a'j.)* to the nearest complex integers 
{{«)*}.=i}Li such that 



i4r 



argmm |aj 



(ar)1,Vz,n. 



(39) 



Remark 5: Qualitative Comparison of the Algorithm with 
Conventional Approaches): 

• (a) Complexity comparisons: The proposed method has 
a complexity that is close to the conventional Distributive 
lA method using alternating optimization [5| but as 
expected converges slower as seen in Table J] This is 
because the subproblems to be solved in each step of the 



proposed algorithm (e.g. Problem (|27t and Problem (|36) ) 
are convex and hence, there exists efficient algorithms 
(such as interior-point method [28]). 

« (b) Limitations of other existing approaches: For the 
conventional lA method |3|, there is always a feasibility 
problem for quasi-static interference channels for K > 
3. All the other baselines (see Section IVIIl i have poor 
performance when K > 2 even when the CSI is perfect. 

• (c) Imperfect CSI considerations: Furthermore, none of 
the baselines consider imperfect CSI while the proposed 
method takes into account imperfect CSI in the design. 
Hence, they have very different performance in the pres- 
ence of imperfect CSI. ■ 

VI. Quantitative Analysis of the Proposed Robust 
Lattice Alignment Solution under Symivietric 
Interference Channels 

In this section, we shall illustrate analytically the potential 
benefits of the proposed robust lattice alignment solution 
versus the brute-force two-stage ML decoding with Gaussian 
random inputs. To obtain a first order comparison, we consider 
a symmetric interference channel as in [1 IJ , where each user 
has one antenna and each transmitter tries to transmit one data 
stream for the desired receiver The symmetric channel model 
is as follows: 



yfe = Xfc ■ 



K 

E - 



■Zfe, 



(40) 



where we assume that the interference channel coefficients h G 
Z + jZ. In im . only real channel coefficients are considered. 

A. Performance of the Proposed Method for Symmetric Inter- 
ference Channels 

Since all channels are symmetric and M = 1, we have 
(vk)* — 1 for all k, and we can consider user 1 without loss of 
generality. During stage I decoding, the optimal interference 
quantization coefficients are given by a* — [0,1,1,...,!]^ 
from Lemma |5] As a result, we would decode the aggregate 
interference [x2 + X3 + ... + x/^J Mod A without any quan- 
tization approximation. By solving the optimization problem 
in ( I27I) . we can obtain the optimal first stage decorrelator, 

^ ^ V PJ 1 + (if- l)|/l|2 + 1/P' 



(41) 



where W = [1, /i, • • • , h]^ , and hence we have 
i?2 = Rs- 



K-1 + {K - 1)P 

During the second stage decoding, we have — h and 
(uj,)* = 1, where the data rate for user 1 is given by Rl = 
log(P). As a result, the minimal achievable data rate is given 
by: 

■ 1 + [{K - l)\h\'' + 1]P\\ 



R„ 



min log(P), log 



K-1 + {K-1)P 



(43) 
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TABLE I 

Comparison of matlab simulation time and the converged performance for a channel realization at transmit SNR=1.5dB for 

m = n = 2 and l = 1 with perfect csi e = 



Number of Users 


Methods 


Time (s) 


Worst-Case Goodput (b/s/Hz) 


Sum Goodput (b/s/Hz) 


K -2. 


Proposed Method 


4.442407 


1.4864 


4.4593 


Distributive lA 


1.579973 


0.3306 


2.4724 


K - 4 


Proposed Method 


9.009197 


0.9036 


3.6144 


Distributive lA 


2.995086 


0.2537 


2.7012 



As shown from (l43T l. the performance bottleneck of the 
proposed method is the stage I decoding at the medium SNR 
regime. Furthermore, the proposed method achieves similar 
performance as in 111] and hence, the proposed method 
is backward compatible with that in ifTTl under the same 
symmetric interference channel realizations. 



B. Comparison with Brute-Force Two Stage ML Decoding 
with Gaussian Inputs 

In this section, we shall compare the performance of the 
proposed method with a baseline, namely the brute-force two- 
stage ML decoding with random Gaussian inputs. Specifically, 
each transmitter transmits random Gaussian signals and 
each receiver of the A'-user symmetric interference channel 
in ( |40l l performs brute-force two-stage ML decoding. With- 
out loss of generality, we consider user 1 in the analysis. 
During the first stage decoding, the Gaussian interferences 
{x2,...,x/f} are decoded by using brute-force ML. The de- 
coded aggregate interference is canceled from the received 
signal yi and the desired signal xi is decoded by a second 
stage ML decoding. Note that unlike the proposed method, 
random Gaussian signals are transmitted and hence, the inter- 
ference space is random. In order to have successful decoding 
in both stages of this system, the data rate has to satisfy the 
following constraints [9J: 



Ri < log (1 + P) 

_ / {K-l)P\h\^ 
E...^^ ^ ^+ 1+V 



(44) 



Therefore, the minimal achievable data rate for all the users 
is given by: 



Rb2 = 

™n(p^log(l 



{K~l)P\h\ 
1 + P 



,log(l + P) , 



(45) 



where the bottleneck is also to decode the interference at 
medium SNR regime. Similarly, the conditions in which the 
proposed method offers a performance gain over the baseline 
system is given by: 



I K-1+{K-1)P 



> 



log 1 



{K-l)\h\-'P 
1 + P 



(46) 



The condition reduces to \h\'^ > (1 + -p) for large K and 

K-l 

> ~^ for medium SNR. In other words, only 

mild conditions are needed for a performance gain over the 
basehne system. 



VII. Performance of Proposed Method under 
AT-usER Interference Channels 

A. Comparison with Conventional Interference Alignment So- 
lution 

As pointed out in Section IIII-AI conventional lA can 
achieve DoF in the 3-user quasi-static MIMO interference 
channels with M > 1 antennas for all transmitters and 
receivers (3). However, it requires full rank channel matrices 
whose coefficients are randomly drawn from a continuous 
distributiorll, as well as the feasibility condition. On the other 
hand, our proposed method works for any channel realizations 
and does not have feasibility problem. In this section, we shall 
illustrate that the proposed method could achieve the same 
DoF as the conventional lA [3] if it is feasible. Without loss 
of generality, consider the case when M is even and denote 
the solution of the alignment method in |3| by 

V,.(IA) = [vi(IA),...,vf/^(IA)], 



U,(IA) = [ui(IA),...,uf/2(IA)] 



(47) 



where v^(IA) and u5.(IA) are the precoder and decorrelator 
for the Z-th data stream of user k respectively, such that 

ui(IA)Hfcfcvi(IA)^0; 



<(IA)HHvf (lA) = 0,Vi 7^ fc or Vn ^ I. 



(48) 



The same DoF of ^ can be achieved in our proposed 
method by setting L = ^ and choosing the design parameters 
as 



4 = 0;4 = 0;4 = l;v^=vt(IA); 

)||2 + p||J|2 u'j,(IA), 



(49) 



1^ l|uL(IA)l 

where ( {u[{lA))"Hkkv{{lA) 0. Since a'^^ = 0, there 
is no stage I decoding and hence, from ( [TtI i. the data rate 
constraint on the data stream is oo. This means that there is 
no data rate constraint associated with stage I decoding. As a 
result, the desired data stream xj. is decoded directly at stage 
II. From JTSl l and ( |49] l. the achievable data rate i?^, for data 
stream xj, is given by: 



Ri < log ( 1 



(50) 



|ui(IA)|P^ 

Therefore, the DoF 3L = can also be achieved in 
our proposed method. Hence, our proposed method achieves 
the same DoF performance as the conventional interference 
alignment method in 1 3 1, |4|, as long as the problem is feasible. 
On the other hand, when the conventional alignment method 
is infeasible, our proposed solution still offers significant 
performance gains compared with various baseline systems. 

'Therefore, the special cases are not considered in conventional lA. e.g., 
key hole effect, line-of-sight and closely spaced anteimas. 
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B. Numerical Simulations 

In this section, we shall compare the system goodput (b/s/Hz 
successfully received) of the proposed robust lattice alignment 
method (designed to optimize the worst-case data rate) with 
five baselines. Baseline 1 refers to the TDMA method. Base- 
Une 2 refers to the brute-force two-stage ML decoding with 
Gaussian inputs described in Section [VI-B| Baseline 3 refers to 
the generalized HK method for K > 2 users, where each user's 
message is split into one private part and several common 
parts as described in [12]. Baseline 4 refers to distributive lA 
based on the alternating optimization method in @, which 
tries to minimize the sum leakage interference. Baseline 5 
refers to the conventional lA method in [3J. In the simulation, 
we assume that all the channel coefficients are i.i.d. zero mean 
unit variance circularly symmetric complex Gaussian. 

Fig. |6] and Fig. [T] illustrate the average worst-case goodput 
and the sum goodput versus transmit SNR (dB), respectively, 
at CSI en-ors e = {0,0.1} for K = 3, M = N ^ 2 
and L = 1. It can be observed that our proposed method 
outperforms all the baselines under both perfect and imperfect 
CSI. Furthermore, the conventional lA and distributive lA 
methods achieve good performance with perfect CSI but they 
are very sensitive to CSI errors, especially in the high SNR 
region. On the other hand, our proposed method is robust to 
CSI errors and outperforms all the baselines as illustrated in 
Fig. |6] Note that our optimization objective is chosen to be 
the worst-case data rate for fairness consideration, which may 
cause some penalt}0 on the system sum rate (as illustrated in 
Fig. Eta)). 

Fig. m illustrates the average worst-case goodput versus 
transmit SNR (dB) with CSI errors e = {0,0.1} for if = 4, 
M ~ N ~ 2 and L — \. Note that in this configuration, the 
conventional lA method |3| is not feasible but the proposed 
method could achieve significant performance gain over all the 
other baselines. 

Fig. |9] illustrates the average sum and worst-case goodput 
versus CSI error e at transmit SNR=11.5dB for K = 3, M = 
N ^ 2 and L = 1. It can be observed that all the baselines are 
very sensitive to imperfect CSI. Even a small CSI error will 
result in a significant degradation in the system performance. 
On the other hand, our method is quite robust to imperfect 
CSI. Fig. [To] illustrates the cumulative distribution function 
(CDF) of the worst-case instantaneous mutual information per 
user at the transmit SNR=11.5dB for K = 4, M = N = i 
and L — 2 with CSI error e — 0.1. Observe that our proposed 
method still offers significant performance gain in the outage 
sense. 

Fig. [TT] illustrates the average worst-case goodput per user 
versus the number of users K at transmit SNR=1.5dB for 
M = N = 4 and L = 2 with CSI eiTor e = 0.1. Observe that 
the performance degrades as K increases, and it is more dif- 
ficult to decode the desired signal. Nevertheless, the proposed 
method outperforms all the baselines at all K. Finally, we 
compare the performance of the proposed method at low SNR 



Had we changed the optimization objective to sum rate, the proposed 
method would also outpeiform the others in Fig. (TJa). 



a 5 

o 

S 3 
< 

1 



Baselines (conventional lA) 
Baseline 4 (distributive lA) 




Baseline 2 
(brute-force two-stage ML) 

_ _ - - - > > 

^ \ 

Baseline 3 (HK method) 



5 1 



5 10 16 20 

Per-User Transmit SNR (dB) 

(a) Perfect CSI (CSI error e = 0) 



Baseiine 4 (distributive lA) 



Proposed metliod 




- - - 
Baseline 1 (TDMA) 



Baseline 3 (HK mettiod) 



10 15 20 

Per-User Transmit SNR (dB) 



(b) Impeifect CSI (CSI error e = 0.1) 



Fig. 6. Comparison of average worst-case goodput (b/s/Hz successfully 
received) versus transmit SNR (dB) with CSI eiTors e = {0, 0.1}. The setup 
is given by K = 3 (number of users), M = N = 2 (transmit and receive 
antennas) and L = I (number of data stream). 



with a naive method that treats interference as noise^ in Fig. 
[T2I Observe that the proposed method outperforms the naive 
method at low SNR and the gain is contributed by precoding. 

VIII. Conclusion 

In this paper, we proposed a robust lattice alignment de- 
sign for K-usei quasi-static MIMO interference channels. 
We exploit the structured interference and propose a robust 
lattice alignment method for irrational interference channels 
with imperfect CSI at all SNR regimes as well as a two- 
stage decoding algorithm to decode the desired signal from 
the structured interference space. For fairness, we maximize 
the worst-case data rate by formulating the design of the pre- 
coders, decorrelators, scaling coefficients and the interference 
quantization coefficients as a mixed integer and continuous 
optimization problem. By using an alternating optimization 
technique and by incorporating imperfect CSI in the design, 
we derive robust solutions with low complexity. Numerical 

'it is shown in 1 2 1 that treating interference as noise is a reasonable strategy 
at low SNR. Furthermore, the naive method does not require CSI at the 
transmitter and is a robust method. 
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Fig. 7. Comparison of the average sum goodput (b/s/Hz successfully 
received) versus transmit SNR (dB) with CSI etTors e = {0, 0.1}. The setup 
is the same as that of Fig. [6] 



Fig. 8. Comparison of average worst-case goodput (b/s/Hz successfully 
received) versus transmit SNR (dB) with CSI eiTors e = {0, 0.1}. The setup 
is given by K = A (number of users), M = N = 2 (transmit and receive 
antennas) and L = 1 (number of data stream). Note that in this configuration, 
the conventional IA method is not feasible (T). 



results verify the advantages of the proposed robust lattice 
alignment method compared with various baseline systems. 



estimate of the Tj, at receiver k is given by: 



Appendix A 
Proof of Lemma[3] 

From (l2Tl i. the real part of the structured aggregate interfer- 
ence is given by: 



(E.,„W}tr-3KH?)Mod A- 
E.,„5i{and?-3{afKlMod A. 



(51) 



Thus, decoding is equivalent to decoding T\, 



E»,„ (5R{<}tr - 9{<}tr 



Mod A. Note that the rate of 



ij, is determined by the maximal rate of Va" 7^ 0. Without 
loss of generality, the maximal data rate is denoted as i?max 
and the corresponding nested lattice is given by A,„ax- The 



Qa.„. (Tj,+z'fc)]Mod A, 



^■{a?}dl 



Mod A 
(52) 



where 



xs{{u[)"Uk^vf - af}3{xf}. From Lemma [1] 5R{xf} and 
3{x"} are independently uniformly distributed over V with 



= P/2, and E[xf] = 0. 



From Lemma|2] the density is upper bounded by the density of 
an i.i.d. zen 
approaches 



an i.i.d. zero mean Gaussian vector zj^, whose variance (ct^,)^ 



P 



^|(ui)^H,,vr 



(53) 

We set the volume of Vmax as Vol(Vmax) > (27re(CT[.)^)^^^, 
since A[, is AWGN good, the probabihty of Pr{T[, ^ T[} 
goes to zero exponentially in T lfT6l - llT8l . Let G(A) denote 
the normalized second moment of lattice A, then Vol(V) — 
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Fig. 9. Comparison of average sum and worst-case goodput (b/s/Hz 
successfully received) versus CSI error e at transmit SNR=11.5dB. The setup 
is given by ii" = 3 (number of users), M = N = 2 (transmit and receive 
antennas) and L = I (number of data stream). 
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Fig. 11. The average worst-case goodput per user versus the number of 
users K at transmit SNR=1.5dB for M = N = 4 and L = 2 with CSI error 
e = 0.1. 
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Fig. 12. Illustration of worst-case goodput (b/s/Hz successfully received) at 
low SNR with CSI errors e = {0, 0.1} for = 4, M = TV = 2 and L = 1. 
The naive method refers to treating interference as noise without requiring 
any CSI at the transmitter. At low SNR, treating interference as noise is a 
reasonable strategy (2] and is also a robust method as no CSIT is needed. 
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Fig. 10. The cumulative distribution function (CDF) of the instantaneous 
worst-case mutual information per user at the transmit SNR=11.5dB for K = 
4, M = 7V = 4andL = 2 with CSI error e = 0.1. 



EH- From ( fT3] l. the rate of the nested lattice code 

max 

IS given by: 



T VVol(V,nax) 



P/2 



G(A)2^e(a[) 

Furthermore, for arbitrarily small ^ > 0, and large enough 
T, G{A)2Tre < (1 + IMl- By Lemma |2l if T is large 
enough, (cr[)^ < (1 +0-^1- ^ result, the following rate is 
achievable for Rmux' 



(54) 



1 



Ids 



P 



|uilP + ^E..JK)^Hfe.vr 



-21og(l+e). 
(55) 



By choosing ^ as sufficiently small, we can neglect the second 
term, and the first term can be approached as desired. Sim- 
ilarly, can also be successfully decoded from 5{y^} 
with the same achievable data rate in ( ISSl l for R^„^. Note 



that, if ai 



5± (all 



except that 



•'max- 



1), 
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= [tf„] Mod A = t^j, and Mod A = tf„. Therefore, 

xf„ can be decoded in stage I and completely nulled out at 
stage II. 

Since the data rate i?" is twice the lattice codes £" (due to 
the real and the imaginary parts), the minimal achievable data 
rate of i?" due to the effect of the CSI error is given by 



minA 



-'P2:,.,j("L)"(Hfc,-A,o 



(56) 
is to choose 

A hi maximizing each term of the sum in the denominator, 
which is given by: 



For given {ajj., u^, v^}, a lower bound of 



K)^(Hfc,-Afc,)vr-ar 



,1 \H 



n|2 



(57) 



where 



< 



< 



lK)^A,,vrP = Tr{(ul.)^A,,v«(vr)«A|^u^ 
Tr|ui(ui)«A,.,vr(vf)^A^} 
Tr{4(ui)^}Tr{Afe,vf(v?)^A^] 
Tr{uiK)^}Tr{vr(v?)^A^A 
Tr{uiK)^}Tr{vr(vf)^}Tr{Af,Afc,} 



J l|2| 



12,2 



Therefore, i?niin given in (d6] i is lower bounded by: 



> 



log ( — 

V I u' 



(58) 



(59) 



which shows that i?" < /i^ is the sufficient condition for the 
achievable data rate Rf under imperfect CSI. 

The proof of the Lemma |3] is an extension of the results 
in l,18J . However, there are the following differences. Firstly, 
the results in ifTsl are for perfect CSI while our results are an 
extension of lITSl for imperfect CSI. Secondly, the results in 
lITSl consider single stream relay systems while our results 
are extended to multi-stream MIMO interference channels. 
Finally, the results in fTS) are analysis based whereas our 
solution is optimization based where we shall optimize the 
precoders v = {{v^j^j^}^]^, first stage decorrelators u = 



{{^k}^=i}f=i^ '■^^ scaling coefficients c 



{{cilLJiiand 

{{^A:}fc=l}i=l- 



the interference quantization coefficients a 
The detailed formulation is illustrated in Section IIV-EI 

Appendix B 

Convergence Proof of Subalgorithms A and B, 
AND Top Algorithm 1 

We first provide the convergence proof of subalgorithm 
A. Note that Given {v, a}, after each iteration in subalgo- 
rithm A (the alternating optimization between u'j. and c'j.), 
the data rate /Ij, in (fTSl ) is increasing. Specifically, we have 
+ > Miju^["i),4(m)) in step 3 of 
subalgorithm A, and we have fi^, (u'j,(m -I- l),c^(m -I- 1)) > 



fore. 



(m)) in step 4 of subalgorithm B. There- 



lA K("^ + 1), + 1))^ ^ (60) 

MiK(m + l),4(m)) >/i^u'jTO),4(m)). 



As a result, subalgorithm A converges. In addition, we have 
formulated subalgorithm B to a standard interior-point algo- 
rithm, i.e., the barrier method [28 j . Therefore the results of 
subalgorithm B {v*, a*} converge to a (^-optimal solution for 
fixed {u,u, c}, and are subject to the tolerance of <, [28 1. 

Note that i?min — niin/c,; (/^fc, ^^d by using the 
above results we have /i'^. (u'j. (m + 1), v(to), a(rn)) > 
MfeKM'VM,a(TO)), and JI^u^(TO + l),4(m + 
1), v(TO),a(TO)) > ^\{u\{m),c\{m),^r{m),a{m)) 
in step 2 of top algorithm 1, in other words, 
-Rmin(u(TO + l),u(m + l),c(m + 1), v(m), a(m)) > 
_Rmin(u(TO), u(m), c(m), v(to), a(m)) . Furthermore, we 
have t(u(rn + 1), u(m + l),c(m + l),v(m + l), a(m + 1)) < 
t{u.{m + 1), u(to + 1), c(m -|- 1), v(m), a(TO)) in step 3 of 
top algorithm 1, where t is given in ( |36] |. In other words, 
-Rmin(u(TO + 1), u(m + 1), c(m -I- 1), v(m + 1), a(m + 1)) > 
i?min(u(m + 1), u(to + 1), c(m + 1), v(m), a(m)) . Therefore 
we can conclude that 

-Rmin(u(m + 1), u(m + 1), c(m + 1), v(m + 1), a(m + 1)) 

> -Rmin(u(m + 1), u(to + 1), c(m + 1), v(to), a(TO)) 

> i?i„in(u(m), u(m), c(m), v(to), a(TO)) . 

(61) 

Since objective i?i„in increases after each iteration as indi- 
cated in (l6Tl i. Algorithm 1 shall converge. 
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